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Kis form factor with Nj = 2 + 1 dynamical domain wall fermions: A progress report 
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We present the latest results from the UKQCD/RBC collaborations for the K13 form factor with 
2 + 1 flavours of dynamical domain wall quarks. Simulations are performed on 16^ X 32 X 16 and 
24"^ X 64 X 16 lattices with three values of the light quark mass, allowing for an extrapolation to 
the chiral limit. After interpolating to zero momentum transfer, we obtain the preliminary result 
/+ "(0) = 0.9609(51) (or A/ = -0.0161(51)), which is in agreement with the result of Leutwyler & 
Roos. 



I. INTRODUCTION 

K — > ttIv (Kis) decays provide an excellent avenue 
for an accurate determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) [ij quark mixing matrix element, \Vus\- 
This is done by observing that the decay amplitude is 
proportional to |T4sP|/+(9^)P, where /+((7^) is the form 
factor defined from the K ^ tt matrix element of the 
weak vector current, = 57^ u 

{7:{p')\V,\K{p)) {p,+p'^)f+iq^) + {p,-p'^)f^{q^) , 

(1) 

where ~ {p — p')^. 

In chiral perturbation theory (ChPT), /+(0) is ex- 
panded in terms of the light pseudoscalar meson masses, 
771^, rriK, m,-i 

/+(0) = l + /2 + /4 + ..., (/« = 0(m:,^^„)) . (2) 

Conservation of the vector current ensures that /+(0) = 1 
in the SU{'i) flavour limit. Additionally, as a result of 
the AdemoUo-Gatto Theorem [2], which states that /2 
receives no contribution from local operators appearing 
in the effective theory, /2 is determined unambiguously in 
terms of mjr, rriK and /^r, and takes the value /2 = —0.023 
at the physical values of the meson masses Q- 
Our task is now reduced to one of finding 



A/ = /+(0)-(l + /2) 



(3) 



In order to obtain a result for /+(0) which is accurate 
to ~1%, it is sufficient to have a 20-30% error on A/. 
Until recently, the standard estimate of A/ = —0.016(8) 
was due to Leutwyler & Roos 3 , however a more recent 
ChPT analysis favours a positive value. A/ = 0.007(12) 
[4] . A calculation of A/ on the lattice is therefore essen- 
tial. 
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The last few years have seen an improvement in the 
accuracy of lattice calculations of this quantity [^, 0, 
H 0] J with the results favouring a negative value for A/ 
in agreement with Leutwyler & Roos. 

The UKQCD and RBC collaborations have embarked 
on a program to perform the first unitary (i.e. Nf — 2-\-l 
fiavour) lattice calculation of the if — > tt form factor us- 
ing dynamical domain wall fermions at light quark masses 
and on large volumes. Preliminary results from this 
study were introduced in Ref. [1^]. The work presented 
here represents an improvement on the earlier analysis 
through a more in depth analysis of the systematic errors 
involved in the extraction of /+(0) from lattice correla- 
tion functions. 

II. LATTICE TECHNIQUES 

A. Parameters 

We simulate with Nf = 2-1-1 dynamical flavours gen- 
erated with the Iwasaki gauge action [ll[ ai (3 — 2.13, 
which corresponds to an inverse lattice spacing « 
1.62(4) GeV [l^, and the domain wall fermion action 
[ist with domain wall height M5 = 1.8 and fifth di- 
mension length Ls = 16. This results in a residual 
mass of amies ~ 0.00308(4) [l2|. The simulated strange 
quark mass, arris = 0.04, is very close to its physical 
value [ij , and we choose three values for the light quark 
masses, aniud = 0.03, 0.02, 0.01, which correspond to 
pion masses « 630, 520, 390 MeV [11]. The cal- 
culations are performed on two volumes, 16^ x 32 and 
24^^ X 64, at each quark mass. See [l^ for further simu- 
lation details. 

B. Extracting form factors 

The techniques we use to calculate /+(0) are similar 
to those originally proposed in Q and have been out- 
lined in detail in Ref. 0]. We restrict ourselves here to 
highlighting the main points. 
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FIG. 1: Ratio for /o(gmax): as defined in Eq. ((Sjl, 

for three simulated light masses arriud ~ 0.03, 0.02, 0.01 on a 
24^ X 64 volume. Further simulation parameters can be found 
in [12. 

We start by rewriting the vector form factors given in 
Eq. ID) to define the scalar form factor 



(4) 



which can be obtained on the lattice at qf^^y. = (m 
m^)^ with high precision from the following ratio [l 



R{t',t) 



Cf"ft^^;0,0)CI^ft^^;0,0) 



where the three-point function is defined as 



{0\OQ{t\y)\Q){Q\V^{t,x)\P){P\oU0)\0) 



(5) 



(6) 



with P, Q ~ IT OT K and Ojr{K) is a local interpolating op- 
erator for a pion(kaon). We note that R{t' , t) — 1 in the 
S'[/(3)fiavour symmetric limit, hence any deviations from 
unity are purely due to SU (3)flavour symmetry breaking 
effects. 

In Fig. [1] we display our results for R{t' , t) for each of 
the simulated quark masses as obtained on the 24'^ x 64 
lattices. 

It is immediately obvious that /o(9max) can be mea- 
sured with a very high level of statistical accuracy. We 
also note that the ratio becomes larger the further we 
move away from the >5'C/(3)fiavour limit. Since there is 
no spatial momentum involved in this ratio, the results 
obtained on the two different volumes should agree, and 
any difference can only be due to finite size effects. We 
find that within statistical errors, finite size effects on 
/o(9max) are neghgible [10]. 

In order to to extract the form factor at = 0, we 
need to obtain results at finite < and by studying 
the dependence we are then in a position to interpolate 




(aq)-^ 

FIG. 2: Scalar form factor fo{q^) for bare quark mass arriud = 
0.02. Results are obtained on two volumes V — 16^ x 32 (red 
diamonds) and V = 24'^ x 64 (blue triangles) . The solid line is 
the result of a fit using a monopole ansatz (Eq. (O) while the 
dotted line is the result of a fit using a zfit ansatz (Eq. (|12p V 



to 5^ — 0. We calculate the form factor /o (9^ < 0) 
using the ratio method outlined in Ref. 0| which has 
the advantage that no knowledge of Zy is required. We 
simulate with two unique choices of the three- momentum 
transfer, q = 2tt/L{1, 0, 0) and q = 27r/L(l, 1, 0), where L 
is the spatial extent of our lattice. We take advantage of 
rotational symmetry to improve the signal, i.e. we use all 
6 and 12 permutations of (1,0,0) and (1,1,0), respectively. 
The ratio method also allows /o(?^) to be extracted from 
the TT — > iiT matrix element which, due to the difference 
between the tt and K masses, doubles the number of q^ 
values. Finally, since (spatial) momentum on the lattice 
is proportional to the inverse (spatial) lattice size, our 
simulations on two different volumes provides access to 
additional g^. 



III. RESULTS 
A. Interpolation to = 

We present our results obtained on each volume for 
fo{q^) in Fig. O for quark mass aniud — 0.02. In the in- 
termediate g^-range we see good agreement between the 
results obtained on the two different volumes, indicating 
that finite size effects are negligible, at least for the quark 
masses considered here. This means that we now have 
results over a large range of q^ to fit to. 

It is not immediately clear what is the best method 
for fitting /o('?^), however previous work 0, Q has found 
that lattice results are described very well by a monopole 
form 



/o(g') 



/o(0) 



(l-q2/A/2) 



(7) 



where M is the monopole mass. From Fig. [2] we see 
that our data is also described well by Eq. [7] (solid line). 
From this fit, we are able to calculate /o(0) for our three 
simulated quark masses. 
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(m^+m2) [GeV^] 
FIG. 3: Chiral extrapolation of A/ using Eq. 0. 



B. Chiral Extrapolation 

Now that we have obtained results for /+(0) = /o(0) at 
three different quark masses, we are in a position to at- 
tempt an extrapolation to the physical pion mass. Insert- 
ing these results into the expression given in Eq. ([3]) , to- 
gether with /2 calculated at the simulated quark masses 
using the ChPT formula [1, HBli we are now left with the 
task of chirally extrapolating A/. 

The Ademollo-Gatto Theorem implies that to leading 
order A/ oc (m^ — rriud)^ , hence we attempt a chiral 
extrapolation using 



A/ 



B{ms - mudY 



(8) 



Note that in the S'[/(3)fiavour limit, A/ = 0, so we expect 
that a fit to our data should produce a « 0, and indeed 
we find this to be the case. In the chiral limit we find 
A/ = -0.0146(28). 

In order to attempt to take into account higher terms 
in the chiral expansion of A/, it has been noted that it is 
convenient to consider an extrapolation of the ratio 0, 



i?A/ 



A/ 



mil 



+ H^K + ml) 



(9) 



We show the extrapolation using this form in Fig [3] from 
which we extract a result at the physical meson masses 
(vertical dotted line). A/ = -0.0161(46). 

C. Estimating Systematic Errors 

As we have seen in the previous section, with only three 
quark masses it is difficult to distinguish between differ- 
ent chiral extrapolation ansatze. Hence, we quote the 
result obtained using Eq. and use the difference be- 
tween the results obtained from the two extrapolations 
(0.0015) as an estimate of the systematic error due to the 
chiral extrapolation. 

Another possible source of systematic error comes from 
our choice of a monopole fit for the dependence of 
/o(<z^). In order to estimate the systematic error due to 
this choice, we fit fo{q'^) at each quark mass with a linear 



form 

/o(9')-/o(0) + alg^ (10) 

a quadratic form 

foiq^) ^ fo{0) + aiq^ + a2q\ (11) 
and a parameterisation proposed in Ref. [l6j 

^ oo 

where |z| < 1 is a mapping oi t ^ q^ onto the unit circle 
in the complex plane 



and t± = (rriK ^rriT^Y . We use io = ^+(1 — \/l — t^/t+) 
as the point that maps onto z = 0, and 



3M- z(<,0) z(t,-Q^) f z{t,to) 



327r 



t 



-t 

-1/4 



to-t 



-1/2 



with = and a^^ — 1.62 GeV is the lattice spac- 
ing. We note that employing different values for to and 
simply affects the fitted values of the coefficients, a^, 
and not the result for /o(0). 

As an example, the dotted line in Fig. [2] shows the 
result of a fit to foiq^) for quark mass arriud = 0.02 us- 
ing the parameterisation in Eq. (|12p truncated to include 
terms up to A: = 2. 

We find that all four parameterisations describe our 
data reasonably well, except the linear ansatz which leads 
to a large /dof on the am — 0.01 dataset. We take the 
difference between the result from the pole fit and the z 
fit as the estimate of the systematic error due to the (f' 
dependence. 

In Fig. 21 we provide a comparison of the value of A/ 
at the physical meson masses as obtained from all four 
parameterisations. The blue triangles indicate the results 
using Eq. ([5]) for the chiral extrapolation, while the red 
diamonds use Eq. ([5]). The four sets of points are, moving 
from left to right, linear (fTO]) . quadratic (|lip . z and 
pole ([7]) fits, respectively. We observe that the systematic 
error due to the q^ dependence is comparable to that 
from the chiral extrapolation, and both are small when 
compared to the statistical error. Hence, one might hope 
that these systematic errors can be reduced along with 
the statistical error by improving the statistics. 

Finally, since we only have results at one lattice spac- 
ing, we are unable to extrapolate to the continuum limit. 
However, lattice artefacts are formally of ©(a^Ag^^); as- 
suming Agp^ ~ 300 MeV we, tentatively, estimate these 



4 



p 
-0.005 - 



-0.02 - 
-0.025 - 




-0.03 I 1 

FIG. 4: Comparison of A/ at the physical meson masses from 
all four parameterisations. Blue triangles indicate the results 
using Eq. ^ for the chiral extrapolation, while the red dia- 
monds use Eq. ((Ojl. The four sets of points are, moving from 
left to right, linear, quadratic, z and pole fits, respectively. 

at w 4%. Hence our preliminary result is 

A/ = -0.0161(46)(15)(16)(7) => /+ ''(0) = 0.9609(51) , 

(13) 

where the first error is statistical, and the second, third 
and fourth are estimates of the systematic errors due to 
the chiral extrapolation, dependence and lattice arite- 
facts, respectively. This result agrees very well with the 
old value of Leutwyler & Roos [3] and earlier lattice cal- 
culations 

While this result for f+{0) is smaller than our earlier 
preliminary analysis (0.9680(16)) fl^ and has a larger 
error, we believe that through a more rigorous investiga- 
tion of the systematics involved in extracting this result, 
we are now much closer to having a finalised result which 
we plan to publish soon. 

Using \Vusf+{0)\= 0.2169(9) from the experimental 
decay amplitude |17| : 

iKsl = 0.2257(9)exp(12)/+(o) (14) 

wefind |K<i|' + |K.|' + |KbP = 1-^, (5 = 0.00076(62), 
compared with the PDG(2006) resuh, 6 = 0.0008(10). 



IV. SUMMARY AND FUTURE WORK 

We have presented a preliminary result for A/ = 
/+(0) — (1 + ./2) using Nf = 2 + 1 dynamical domain wall 
fermions with three choices for the light quark masses. 
Our result A/ = -0.0161(46)(15)(16)(7) agrees very weh 
with the Leutwyler & Roos result Jj and confirms the 
trend of other lattice results BlSi, H which prefer 
a negative value for A/. We performed our simulations 
with matched parameters on two volumes and we observe 
no obvious finite size effects. 

We are in the process of improving this result by 
looking at ways of decreasing the error on the point at 
aiTiud = 0.01. We will also soon have a result at a lighter 
quark mass {aniud = 0.005) which will assist in improv- 
ing the chiral extrapolation. Additionally, this result has 
been obtained at a single value of the lattice spacing, so 
future simulations will need to be performed at least at 
one more lattice spacing to investigate scaling behaviour. 
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